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ABSTRACT 



The Halo Galactic Globular Clusters (GCs) lying inside the solar circle {Inner , Rqc < SKpc, 
[Fe/H] < —0.8) are shown to keep no record of the strong correlation between GC absolute 
integrated magnitudes and core parameters (concentration and central density) clearly de- 
tected for the Halo GCs lying outside this circle {Outer ) as well as (at lower statistical 
significance) for those belonging to the Disk population {[Fe/H] > —0.8). 

We present here both observational evidences and theoretical simulations supporting 
the hypothesis of a primordial origin for such correlations and of an environment-induced 
destruction of them if absent. 

In particular, this hypothesis lead us to argue that Inner Halo GCs have undergone a 
chaotic, environment-driven, dynamical evolution not shared by the Disk GCs though the 
radial zone of the Galaxy occupied by the two groups is nearly the same. 

1. INTRODUCTION 

The dynamical evolution of globular clusters (GCs) is likely to be heavily affected by the 
environment of the host galaxy. In the Galactic GC system, environmental effects are due 
to interactions between the GCs and the Galaxy, either via evaporation mostly controlled 
by the Galactic tidal force at the cluster perigalacticon (Aguilar 1993) or via the dynamical 
shocks the GCs undergo during their passages through (or near) the main substructures of the 
Galaxy, i.e. the disk, the bulge or big gas clouds. All these mechanisms as well as dynamical 
friction can lead to GC destruction (see Aguilar, 1993 for an extensive discussion). 

However, the existing globular clusters have survived nearly for an Hubble time, and the 
expected present-day destruction rate is 5±3 cluster per Gyr (Hut & Djorgovski 1992). This 
implies that only ~ 3% of the whole Galactic GC system will disappear in the next Gyr or 
so. As a consequence, the bulk of the existing Galactic globulars consists of stellar systems 
which have continuously evolved, for several core-relaxation times, toward higher and higher 
concentrations and, in some cases, till core collapse (Djorgovski & Meylan, 1994, hereafter 
DM94). 
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The time scales of such an evolution can be very different from cluster to cluster and 
strictly depend on the efficiency of the various mechanisms involved during the whole cluster 
lifetime. The dynamical status of a globular cluster must therefore store information not only 
about its "intrinsic" evolution but also about the effects induced by the environment through 
which it has passed along its orbits. A straightforward example of this kind of information is 
for instance supplied by the evidence that all of the ~30 clusters which are thought to have 
undergone core collapse lie inside the inner 10 Kpc from the Galactic Center (Chernoff & 
Djorgovski, 1989): this is indeed the proof that environmental effects accelerate dynamical 
evolution of the globulars. 

Djorgovski and Meylan (DM94) have presented the most relevant correlations between 
position in the Galaxy and structural parameters for the whole Galactic GC system, show- 
ing also clear hints on the existence of the quoted "memories" of the past GC history. In 
particular, we recall that: 

1) GC structural parameters ( Tc -the core radius, po -the central density, C -the con- 
centration) are correlated between each other and with the cluster integrated absolute 
magnitude My, in the sense that central concentration increases with increasing cluster 
luminosity (see also van den Bergh 1994). 

2) GC structural parameters and positions in the Galaxy display a good correlation, in the 
sense that clusters far from the Galactic Center tend to be looser. 

3) My and Rqc : where Rqc is the GC distance from the Galactic Center, are found to 
be nearly uncorrelated between each other (r = 0.11), if one excludes from the sample 
the five most distant clusters whose origin is uncertain (see discussion below). 

The basic aims of the present work are: 

a) to detect different dynamical evolutionary conditions or trends in (properly selected) 
different groups of Galactic globulars; 

b) to evaluate the relative importance of primordial and evolutionary mechanisms in each 
group. 
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The simple idea at the basis of this analysis is that the correlation between concentration and 
integrated magnitude is mostly driven by internal conditions, whilst the correlation between 
GC structural parameters and, for example, Rgc or Zqp - where Zqp is the height over 
the Galactic plane - must be related to environmental effects on the GC dynamical evolution. 

Within this approach, Rqc and My are considered to be primary parameters, that is 
directly related to quantities settled at the very origin of the globulars {i.e. apogalactica and 
cluster masses, respectively). Moreover, they are supposed to undergo little changing and 
evolution throughout the life of a typical cluster (Murray & Lin 1992). Finally, since we are 
looking for evolutionary characteristics, our analysis is focused on GC core parameters, which 
certainly are the most sensitive to any dynamical evolution (Murray & Lin 1992, DM94). 

2. DATA-SET AND DEFINITIONS 

The data-sets adopted are those listed by Djorgovski (1993) and Trager et al. (1993). Note 
that the absolute values of integrated magnitudes ( |Mv| ) are used throughout the paper 
as they rank clusters nearly as masses do. From the original entries (143), we have actually 
considered only 113 objects as for the other GCs important data were missing or they were 
classified by Djorgovski (1993) to be of low quality. Furthermore, because of their badly 
defined core parameters, we were forced to exclude from our sample the GCs classified "c" 
and "c?" by Djorgovski (1993), candidate members of the group of the so-called Post Core- 
Collapse Clusters (PCC), (see, e.g.. Lugger et al. 1995 for a recent investigation on the core 
sizes of PCC clusters leading to values significantly different from those provided by Trager 
et al. 1993). The data concerning PCCs were used only in the computation of average values 
for [Fe/H] , |My| and \Zgp\/ Rgc (see Tab. 1) and in evaluating correlations between 
parameters other than po : or C , as for example the above quoted one between Rgc 
and \Mv\ ■ 

The final total sample has then been divided into four groups: 
1. Disk Clusters: having [Fe/H] > -0.8 (Zinn 1985). 
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2. Inner Halo Clusters: having [Fe/H] < —0.8 and lying inside the solar circle, adopted 
to be at Rgc — SKpc. 

3. Outer Halo Clusters: having [Fe/H] < —0.8 and lying outside the solar circle, but inside 
4:0Kpc from the Galactic Center. 

4. Extreme Halo Clusters: situated at distances larger than 60Kpc from the Galactic Center 
(as known, no cluster is found between ~ 40 and ~ QOKpc from the Galactic Center, 
Zinn, 1985; Armandroff, 1993). 

The main characteristics of each group are presented in Table 1. While the reasons to 
consider Disk and Extreme clusters as singular families are straightforward (Zinn 1985, 1993; 
Armandroff 1989), the partition of the halo clusters in two subgroups at a fixed Rqc can 
seem rather arbitrary. The aim was that of examining the dynamical status of two subgroups 
which certainly spend most of their life in different environmental conditions though having 
similar (halo-like) kinematical properties. Such a claim can be supported for instance by 
the results of Chernoff et al. (1986) who found that the impact of disk shocking on the 
GC system is highly effective in the radial range 3 — 8Kpc and also by those presented by 
Weinberg (1994) who, on a quite different basis, suggests that "shock heating plays a defining 
role in the evolution of clusters inside the solar circle" . Furthermore, setting the separation 
limit at 8Kpc, we obtain a sample of halo globulars (the Inner one) actually sharing the same 
radial zone of the Galaxy as the Disk ones (with the only exception of the disk cluster Pal 8 
at Rgc — 20.9Kpc), allowing thus a direct "local" comparison. Finally, it has to be stressed 
that we have verified that the results we present here are essentially unchanged if the limit 
is moved 1 — 2Kpc inward or outward. 

3. OBSERVATIONAL EVIDENCES 

In order to decide whether a correlation between two physical quantities is present or absent in 
a specified sample we use the statistically robust non-parametric Spearman rank correlation 
coefficient s (see Press et al. 1992). In calculating this coefficient the absolute values of 
two given parameters (whose distribution laws are often unknown) are substituted by the 
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ranks of their values. The main advantage is that the distribution function of the ranks (i.e. 
1,2,3,4...) is perfectely known (uniform distribution) and one can safely test the significance 
of any deviation of the measured s from an expected value (say 0, i.e. no correlation). 
In particular, given two vectors a and b having dimension N, the parameter 



is distributed approximately as a Student's distribution with N — 2 degrees of freedom, 
indipendently of the original distribution of a and b, and tests the null hypothesis (Hq) that 
the two vectors are uncorrelated (Press et al. 1992). 

An analogous test of significance can be performed also for the ordinary linear correlation 
coefficient -the Pearson's r (Taylor 1984, Press et al. 1992)- but it is valid only if the 
distributions of a and b jointly form a binormal Gaussian distribution around their mean 
values (Press et al. 1992). 

For clarity and completeness, we report in any plot the values of r (which can at lest 
be regarded as an indication), s (which measures the strength of a correlation), and Xg (i.e. 
the percent probability that the two considered vectors are uncorrelated). In particular, the 
quantity 1—Xs can be regarded as the quantitative confidence level of a given correlation. Note 
that the use of ranks instead of absolute values prevents from the any possible overweighting 
of outliers data which can significantly affect the ordinary linear correlation coefl&cient r. 
In conclusion, the presence/absence of any correlations in the data (and their strength) is 
evaluated in a safe non-parametric way. In particular, we will consider a correlation significant 
when Xs < 5% and very significant when Xs < 1%. 

The correlations detected for the GC groups {Disk , Inner and Outer ) in the planes formed 
by the considered structural parameters versus logRoc and \Mv\ are presented in Table 
2. For each sample the measured s and Xa values are reported for each couple of considered 
parameters. 

3.1. Trends with Rgc 




6 



The main indication emerging from inspecting the first three columns of Table 2 is that the 
dynamical status of the Outer clusters, in spite of the large radial range spanned by the 

sample (8 4QKpc), is practically unaffected by their location within the Galaxy, while 

both the Disk and the Inner clusters show some correlation. To give a view of the different 
behaviours of the various sample, we report in Fig. 1 the distribution of the data points in 
the plane logRoc ~ ^ogpo for the Outer , Inner and Disk samples, respectively. The full 
lines represent the linear regression fits to the data. 

The difference between the Outer group and the others is also testified to by the fact that 
the Outer sample as a whole displays a significantly different distribution of concentrations 
than the Inner and Disk ones, centered on lower po ^ind C values. Moreover, this global 
evidence is consistent with the result of Chernoff & Djorgovski (1989) who showed that the 
fraction of PCC clusters increases with decreasing Rqc j and it gives further support to the 
idea that the evolution of GCs is strongly affected by the Galactic environment, being faster 
for clusters close to the Galactic Center where the tidal field of the Galaxy is much stronger. 

3.2. Trends with \Mv\ 

Turning to Table 2, columns 4-6, it is evident that the Outer clusters show a clear-cut and 
very significant correlation between \Mv\ and structural parameters in any of the examined 
planes. Also the Disk clusters display good correlations, though less evident and significant. 
On the contrary, in spite of the strict similarity in the radial range actually spanned with the 
Disk sample, the Inner clusters do not show any statistically significant correlation. Hence, 
we can point out here the existence of a striking difference between the two subgroups of halo 
Galactic GCs. 

In this respect, it is particularly interesting to consider the plane |My| - logrc , where both 
quantities are direct observables, and where the Inner clusters display (if anything) a very 
weak correlation just in the opposite sense than that of the other two groups {i.e. there is 
no possible doubt on the absence of this correlation in the Inner sample). It is also worth 
recalling that while the correlation is clearly present in the Outer sample and possibly in the 
Disk one, the analysis of the global sample (DM94, van den Bergh 1994) shows only a weak 
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correlation in this specific plane, clearly masking the differences associated to the various 
sub-groups. 

In order to present the results in more detail and to demonstrate their statistical sig- 
nificance, we show in Figure 2 the actual distributions of the clusters included in the three 
samples in the \Mv\ - logpo plane. The scale is the same in any panel and covers the whole 
range spanned by the considered Galactic GCs. The full lines are linear regression fits to the 

and r are the statistical parameters described above. 
By inspecting the two upper panels (concerning the Outer and Inner clusters), one could 
argue that the absence of any clear-cut correlation in the Inner sample is simply due to the 
narrower range spanned in both \Mv\ and logpo than the Outer sample. However, as shown 
in the right-bottom panel, this is not the case. In fact, by extracting from the Outer group 
a subsample (Outer -B) sharing exactly the same parameter ranges as spanned by the Inner 
globulars one gets a still highly significant correlation (> 99% confidence level). Finally, it is 
worth noting that the same kind of analysis applied to the \Mv\ - C and \Mv\ - rc planes 
would yield fully consistent results. 

In conclusion, independently of any apparent bias, it seems quite evident that within 
the Outer sample it is possible to detect a clear-cut trend which is not present in the Inner 
one. As discussed below, this fact may have important implications in the understanding of 
GC formation and dynamical evolution. 

3.3. A bivariate analysis 

Before concluding this Section it may be worth reporting the results obtained by applying a 
bivariate analysis to the three groups to find possible linear combinations of the type: 

Q = a{\Mv\ + blogRoc) 

which could optimize significantly the correlation coefficient with respect to the best mono- 
variate analysis, and where Q can be in turn logrc , logpo , C , respectively. The results so 
obtained are different for the three groups. 
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In the Outer sample, no significant difference has been detected with respect to the 
monovariate framework, being each considered plane actually dominated by the strong rank- 
ing induced by \Mv\ ■ 

The Inner group is slightly affected, in the sense that the contemporary consideration 
of Rgc and \Mv\ marginally improves the correlations. 

The Disk sample shows a clear-cut improvement in the correlations, evidence for an 
intrinsic bi-parametric dependence. As can be seen in Figure 3 where the Disk clusters are 
plotted in the plane logpo versus \Mv\ — S.blogRGC, s increases by ~ 0.2 with respect to 
the monovariate case (the plane \Mv\ - logpo ), and the significance level is very high. In 
particular, the relationship is quite remarkably narrowed, and this is a further confirmation 
that Disk clusters recorded the "signals" impressed by both \Mv\ and Rgc ■ 

3.4- The Extreme clusters 

The poorness of the Extreme sample (6 clusters) prevents any significant statistical analysis 
of the kind performed for the other considered groups. Because of their position and chemical 
properties many authors have argued for an independent origin of these clusters with respect 
to the Galactic GC system (Harris 1976, van den Bergh 1983, Zinn 1985). Fusi Pecci et al. 
(1995) have shown that four of them (Pal 14, Pal 3, Pal 4 and Eridanus), which span a 
narrow metallicity range, have plunging orbits and lie on a plane in the sky passing through 
the Galactic Center (see also Majewski 1994). 

Moreover, while no significant difference is found between the luminosity functions of the 
Disk , Outer and Inner samples, a KS-test shows that the probability that the Extreme and 
the Inner -\- Outer samples are drawn from the same parent population in \Mv\ is less than 
0.3%. Note that if the same KS-test is performed comparing just the Extreme and Outer 
samples the probability increases only up to ~ 0.6%, while for example for the Inner and 
Ottter samples the probability is ~ 81%. Hence it seems reasonable to isolate this small group 
of distant globulars and consider them separately. 

For what concerns the specific analysis carried out above, 5 of these Extreme clusters 
lie on the same relations as found for the Outer globulars in the planes "luminosity vs. core 
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concentration", while only one, NGC 2419, greatly deviates, being NGC 2419 very luminous 
and sparse. This fact is not totally at odds with the general trend as NGC 2419, very bright 
and loose, is anyway much denser than the other, less massive. Extreme clusters. 



The observational evidences presented in Sect. 3 ask for interpretation and modeling. In 
this section we describe the main results of a theoretical investigation on the evolution of 
GC systems (see also Vesperini 1994, 1995a,b for further details) which may add information 
useful to properly interpret the presence or lack of correlation between the various quantities 
in the considered GC samples. In particular, the main goal of the present simulations is that 
of determining the role of primordial conditions and evolutionary processes in establishing 
the observed correlations. Due to the uncertainties on the differences between halo and disk 
clusters, particularly as far as age and interaction with the disk are concerned, we have limited 
our analysis to a population of halo clusters. 

The model assumptions and parameters 

In order to follow the evolution of GCs with a large set of different initial values of mass 
and concentration, and located at different distances from the Galactic Center, we have 
used a method introduced by King (1966) and later on applied by Prata (1971a,b), Chernoff 
et al. (1986), Chernoff & Shapiro (1987), based on the assumption that GC evolution can be 
described, at least until the onset of gravothermal instability, as a sequence of King models 
with time evolving concentration, mass and radius. Here, we give just a brief description of 
the method (for further details see the above references). 

Given the total energy Efot-, the total mass M, and the truncation radius r^, a King- 
model describing an individual GC is univocally determined. We shall assume that the 
truncation radius is determined by the tidal field of the Galaxy (see, e.g., Spitzer 1987) and 
defined as 



4. THEORETICAL SIMULATIONS 




(1) 
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where Rqc is the radius of the orbit of the GC around the host galaxy and Mg is the mass 
of the galaxy contained inside that radius. For simplicity, the cluster orbits are taken to be 
circular and lying on planes perpendicular to the plane of the galactic disk. 

After a given time interval At, the change of mass and energy produced by evolutionary 
processes can be calculated and a new tidal radius can then be determined by the tidal 
limiting condition (eq. 1). The three new values of mass, energy and radius identify a new 
King-model with different concentration and different scale parameters describing the cluster 
at time t + At. The total change in mass and energy is given by the sum of the variations of 
these quantities due to all the processes taken into consideration in the study. 

In our investigation we have included the effects on the mass and the energy of the 
cluster due to disk shocking and internal relaxation (see Chernoff et al. 1986 and Vesperini 
1994, 1995a,b for further details). As for the model of the disk of the Milky Way, we have 
adopted the same one used in Chernoff et al. (1986) obtained by fitting with a two-component 
isothermal model the acceleration along the z direction in the solar neighbourhood determined 
by Bahcall (1984) 



withi?o = 8Kpc, Kq = 3.47xl0-^cm s'^, zq = 175 pc, Ki = 3.96xl0-^cm s-^, zi = 550 pc. 
Acceleration varies with the Galactocentric distance according to the surface density that, 
according to the Bahcall-Schmidt-Soneira (1982) model for the Galaxy we will adopt in our 
investigation, falls off exponentially with an exponential scalelength h = 3.5Kpc. Clusters 
are assumed to rotate with a constant velocity V = 210km corresponding to a constant 
ratio Mg/Rgc = lO^Mo/pc. 

The evolution of each cluster is followed until one of these three conditions is satisfied: 
1. Wo > 7.4 -limit for the onset of gravothermal catastrophe (Wq is the dimensionless 
central potential). Both analytical studies and numerical integrations of the Fokker- 
Planck equation (Katz 1980, Wiyanto et al. 1985) indicate that King models with Wq 
larger than this value are unstable against gravothermal catastrophe. Gravothermal 




(2) 
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collapse is driven by a thermodynamic instability and can not be described by the above 
method. 

2. Wq < 0.05 -under this value the system is conventionally considered to be dissolved 
(Chernoff & Shapiro 1987). 

3. t = 1.5 X 10-'^° years (Hubble time). The present era has been reached. 

Our final sample of globular clusters include only those clusters that can still be represented 
as King models after one Hubble time. 

The main limitations of our model are the impossibility of following the evolution of 
clusters after the onset of gravothermal catastrophe and the fact that we consider only clusters 
on circular orbits. In a future work we will extend our analysis by including clusters on 
eccentric orbits but it must be noted that this leads to a significant increase of the parameter 
space of initial conditions. Restriction to circular orbits implies that some caution is necessary 
in comparing theoretical results with observational data particularly for what concerns trends 
with galactocentric distance, as no mixing of orbits is present in our model. On the other 
hand we point out that a detailed comparison of the theoretical results from simulations 
with observational data is beyond the scope of our analysis; the main goal is rather that 
of establishing whether evolutionary processes can modify given initial conditions for a GC 
system toward the observed properties or it is necessary to apply to initial conditions to 
explain them. 

4-2 Evolution of a globular cluster system 

In order to start our simulation to follow the evolution of a GC system we need to set the 
initial conditions: the distribution of orbital parameters, the initial mass function of the GC 
system, the distribution of initial concentrations of the GCs, and the spatial distribution of 
the GCs in the Galaxy. 

The present knowledge of the real initial conditions of the Galactic GC system is rather 
poor, and it is not clear to what extent the present observational properties of the GCs keep 
memory of the initial ones. Thus we have investigated different initial conditions and here 
we describe the results of two of the simulations we have carried out, pointing out the main 

12 



points of interest connected with the issues addressed in the paper (see Vesperini 1994, 1995b 
for the results concerning other properties of globular cluster systems). 

The results described below have been obtained starting from a sample of 1000 clusters 
located on circular orbits with Galactocentric distances between 1 and 20 Kpc, and with 
a spatial distribution in the Galaxy such that the number of GCs per cubic kiloparsec is 
N{Rgc) oc Rqc^ according to what observed for Rqc > ^Kpc (see, e.g., DM94). The initial 
mass function has been chosen to be a truncated power-law, /(M) oc M~^, with lower cut-off 
Miow = lO^-^M© and upper cut-off Mup = 10^ M©. As for the distribution of concentrations, 
C, we have investigated two possibilities: the first ( run I ) with C and logM initially 
correlated 

C = -2.0 + 0.6 log M (3) 

chosen on the basis of the properties of the more massive and distant clusters (but it is 
important to point out that the qualitative conclusions of the theoretical analysis do not 
strongly depend on the exact form of the initial relation between C and logM). The second 
( run II ) starting with C and log M uncorrelated in order to test the hypothesis that the 
correlation between concentration and mass is due to evolutionary processes (Djorgovski 1991; 
DM94). In run II, the initial distribution of C is uniform in the range [0.2 — 1.5]. 

Figure 4a,b show the final C — logM plane for the two simulations actually carried out. 
run I : starting with C — log M correlated. The correlation is essentially preserved even though 
with a large spread (see Fig. 4a) due to the effects of evolution causing the concentration 
to decrease or increase depending on the initial condition of the system; evolutionary 
processes modify the slope of the relation and give rise to the observed scatter. 
run II : starting with C — logM uncorrelated. The final sample does not show any trend 
between the two quantities (see Fig. 4b), confirming the skepticism about the possible 
role of evolutionary processes in giving rise to the observed correlation. 
Figures 5a,b show the C — log M plane for the final sample resulting from the run I for inner 
{Rgc < SKpc) and outer {Roc > SKpc) clusters: it is evident that the initial correlation be- 
tween concentration and mass of the clusters is more efficiently destroyed in the Inner regions 
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where evolutionary processes are more efficient, while for clusters in the Outer regions the 
correlation is stronger since the memory of the initial conditions is preserved. The Pearson's, 
r, and the Spearman's, s, correlation coefficients for the final C — logM relationship for the 
two samples have the same trend present in observational data being r = 0.95, s = 0.94 for 
the Outer clusters and r = 0.57, s = 0.56 for the Inner clusters. 

Also for the C — log Rqc relationship the final result of the simulation {run I) shows 
the same trend present in the observational data. Initial conditions are such that there is 
no correlation between concentration and galactocentric distance neither for inner nor for 
outer clusters. In the final sample, inner clusters display a trend between C and log Rgc in 
the sense of clusters having higher concentrations at smaller galactcocentric distances, with 
correlation coefficients (r = —0.60, s = —0.64). No correlation between C and log Rgc is 
present (r = 0.17, s = 0.24) for outer clusters also at the end of the simulation. The trend 
for inner clusters to have higher concentrations as the galactocentric distance decreases is 
completely originated by evolutionary processes partly as a consequence of the environmental 
effects accelerating the evolution of clusters located in the inner regions of the Galaxy and 
partly as the result of a selective depletion of low-concentration clusters. 
Trends similar to the observed ones are obtained also for the other relationships concerning 
core parameters-mass and core parameters-galactocentric distance. 

Before closing, it is important to point out that the good qualitative agreement between 
theoretical and observational results has been obtained in spite of our conservative assumption 
on the range of Galactocentric distances for Outer clusters in the theoretical sample [8 — 
20]Kpc compared with the observational one [8 — 40]Kpc. By including larger values of Rgc, 
the fraction of clusters in the sample of the Outer GCs preserving memory of initial conditions 
would increase, making the differences between Inner and Outer globulars even larger. 

5. DISCUSSION 

5.1. A "primordial" origin for the correlation \Mv\ -core parameters 
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The new significant result presented here is that the relation standing between absolute 
luminosity (i.e. total mass) and core parameters {i.e. concentration and central density) of 
Galactic GCs appears to be mostly evident in the group of clusters lying in an external region 
of the halo {Outer clusters), whilst it is substantially absent in the inner halo objects {Inner 
clusters). 

The very high statistical significance of these observed relations suggests that, in spite 
of the fact that the members of Outer group seem to have preferentially plunging orbits (van 
den Bergli 1994), they were poorly disturbed by evolutionary dynamical processes. Moreover, 
these clusters do not show any ranking in concentration induced by their distance from the 
Galactic Center, and this weakens the possibility that the detected trend could have an 
environmental origin. 

Taking into account the above evidences and the results of our theoretical simulations, con- 
cerning the origin of the detected correlation luminosity- concentration we are led to conclude 
that: 

The correlation has a primeval origin, in the sense that more massive clusters were born 
more centrally concentrated. In our view, evolutionary processes are likely to have the 
efi"ect of destroying the correlation established at the time of cluster formation. In fact, 
both the observational data and the results of our theoretical investigation show that 
the studied correlations are much weaker for clusters in the inner regions of the Galaxy, 
where all evolutionary processes are more efl&cient. In this respect, it is interesting to 
note that a similar correlation between concentration and luminosity holds for dwarf 
elliptical galaxies (see, e.g.. Fig. 10b in Binggeli et al. 1984, Fig. 11 in Ichikawa et al. 
1986) for which evolutionary processes are likely to be irrelevant. The confirmation 
of this working hypothesis would obviously be of great relevance in constraining the 
physical processes playing a role in the formation of these stellar systems. 
Alternatively, one can hypothize that there is an internal mechanism accelerating the evolu- 
tion of the more massive clusters and making them more and more centrally concentrated. 
Such a mechanism has necessarily to be internal because the mass of the cluster is (at first 
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order) determined only by the primordial conditions at the time of its formation. As a con- 
sequence, this kind of evolution has to affect also an ideal isolated cluster, i.e. in the quoted 
range of parameters a larger mass in some way accelerates two-body relaxation and, in turn, 
the evolution of the cluster toward higher concentrations. Still, for a given size of the system, 
two-body relaxation time is proportional to M-^/^, and this runs just in the opposite direction 
needed to explain the observed trend via evolutionary effects. 

5.2 Comments against the "environmental" origin 

The sample Galactic globulars here considered is necessarily made just of survived GCs. One 
might thus imagine that any relation involving cluster structural properties could well be 
originated by selective depletion processes within the original GC system. 
However, while it is acceptable that the loosest massive clusters were somehow disrupted, it 
seems very difficult to understand which process could have destroyed the more massive tail 
of the "loose cluster" distribution, leaving the lighter population still alive. 

In order to explain the global C - \Mv\ and logpo - \Mv\ trends, DM94 have invoked "... 
a differential survival effect, with more massive clusters surviving longer and reaching more 
evolved dynamical states". We are however skeptical about this scenario mainly because, 
since present-day clusters have all survived till now, the only way to actually allow massive 
clusters to evolve longer is to postulate that less massive GCs are significantly younger than 
the others, which is so far not confirmed by any observational evidence. Moreover, it is 
unlikely that the high-mass, low-C GCs in the outer regions of the Galaxy could evolve 
enough to depopulate completely the high-mass, low-C part of the logM — C plane. As 
shown also in the results of our theoretical simulations ( run II), if present in the primordial 
GC system, high-mass, low-C clusters should still exist at present. 

Finally, assuming that the effect of disk shocks could be important also at large Galac- 
tocentric distances, one could imagine that the GCs more massive than a certain critical 
value (M*) are induced by the shock events to contract their cores, while those less massive 
than M* are induced to expand. This is exactly the way in which disk shocks are believed to 
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act, but the critical parameter is found to be the concentration and not the mass (Chernoff 
et al. 1986; Weinberg 1994). So this is a viable explanation for the detected trend only if a 
7?nmor(im/ "mass-concentration" coupling is postulated. 

5.3 Differential impact of disk shocks in the inner regions of the Galaxy 

Another observational evidence seems to yield a further important indication as the impact of 
environmental effects are apparently different on different groups of clusters, even if located 
within the same region of the Galaxy. 

In fact, the Disk clusters, which actually occupy the same radial zone of the Galaxy as the 
Inner globulars, keep record of the internal processes (being ranked by the quoted linear 
combination of \Mv\ and logRac )i while the Inner clusters do not. What is the cause of 
such a different behaviour for the two groups? 

In our view, the most viable explanation is that, due to their very different orbits, the 
Inner globulars experience stronger interactions with respect to Disk clusters, during the 
passage through the same Galactic environment. Disk clusters are orbiting circularly near 
the plane of the Galaxy and they are rotating with approximately the same velocity of their 
stellar neighbours (Zinn 1985, Armandroff 1989, Majewski 1993), so minimizing the shocking 
effect of the Galactic Disk (Aguilar, 1993). On the contrary, the Inner clusters are known 
to have a small net rotation compared to that of the Galactic Disk (Armandroff 1989, Zinn 
1993), and they have a distribution of the heights on the Galactic Plane spread out over a very 
wide range and suggestive of orbits incident on the Galactic Plane with large inclinations. 
Due also to the rapid Galactic crossing time (Chernoff & Weinberg 1990), they had probably 
undergone a great number of passages through the densest parts of the Galactic disk, suffering 
several episodes of accelerated evolution (DM94) which erased any memory of the ranking 
settled by the primordial formation processes. 

To get a direct confirmation of the different orbital characteristics actually displayed 
by the two groups, we have performed a straightforward test of flattening comparing the 
distributions of the parameter {\Zgp\/ Rgc) in the two {Disk and Inner ) groups. 
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A KS-test rejects the hypotesis that the two samples are drawn from the same population at 
99.9%, the Disk sample being significantly flatter on the Galactic Plane than the Inner one. 
Moreover, looking at the data presented in Table 1, it can also be noted that the Disk sample 
is the only one showing a < \Zgp\/Rgc > value hardly compatible with a uniform spherical 
distribution around the Galactic Center. 

Clues on a possible major effect of disk shocking on Galactic globulars have been recently 
pointed in a previous paper (Bellazzini et al. 1995), where we showed (based on the observed 
frequency of Low Mass X-ray Binaries in the Galactic GCs) that the increasing of the slope 
of GC IMF with increasing Rqc and Zap (Djorgovski, Piotto and Capaccioli, 1993), if 
real, cannot be primordial but it can be settled by the selective depletion of the lighter stars 
from the halo of the clusters during their passage through the densest internal parts of the 
Galactic Disk, as already envisaged by the same latter authors. 

6. CONCLUSIONS 

The results presented here add new observational evidences on the impact of disk shocks 
on Galactic GCs. Further, among the various correlations detected in different GC groups 
between GC integrated absolute magnitudes and intrinsic structural parameters, it is now 
possible to discriminate those likely to be primordial and those likely to result from evolution. 
A proper set of theoretical simulations on the evolution of a GC system reproduce the observed 
trends and support our claims. 

The main conclusions drawn from the present analysis are: 

1) Considering the whole set of Galactic GCs, the already known evidence that GC concen- 
tration increases with increasing absolute luminosity is confirmed. Our analysis shows 
in addition that the correlation is very strong for clusters located in the outer regions 
of the Galaxy, where the tidal field is weaker, and almost absent (probably because de- 
stroyed by evolution) in the inner halo. This lends strong support to the hypothesis of 
a primeval origin of the observed correlation, where it exists. In this respect, the results 
of theoretical simulations on the evolution of GC systems we have carried out show that 
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evolution tends to erase the primordial concentration -total luminosity correlation and 
not to create it as claimed by some authors (Djorgovski 1991, DM94). 

2) The Inner clusters show a much better correlation between core parameters and Galacto- 
centric distance, Rgc , than the Outer globulars. This is likely the result of evolutionary 
processes, more efficient as Rgc decreases. The Outer clusters are probably closer to 
their primeval status, and do not show any significant trend with Rqc • 

3) The impact of the environment has been apparently very different for two groups of 
clusters (Inner and Disk ) actually sharing the same radial zone of the Galaxy. In 
fact. Disk clusters display a well defined correlation, whilst Inner globulars do not. A 
plausible mechanism able to cause such a difference is disk shocking which had a major 
impact on the inclined, high Zqp orbit of the Inner clusters and a neglegible one on 
the very fiat, rapidly rotating Disk system. 

4) The galactic environment plays apparently a major role in the dynamical evolution of 
GCs only in the inner 8 — lOKpc from the Galactic Center, while the population presently 
found outside this circle seems to be almost unaffected by environment-induced effects, 
preserving very well their primeval intrinsic properties. 
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Figure Captions: 

Figure 1. Distributions of tlie tliree cluster samples {Outer, Inner and Disk) in the plane 
logRoc - logpo ■ The statistical parmeters s, Xg and r are described in Sect. 3. The full 
lines represent the best fit linear regression loci. 

Figure 2. Distributions of the three cluster samples {Outer, Inner and Disk) in the plane 
\Mv\ - logpo . The statistical parameters s, Xg and r are described in Sect. 3. The full lines 
represent the best fit linear regression loci. The low-left panel shows the subsample of Outer 
clusters spanning the same range of the observable plane as the Inner globulars. As can be 
seen, the correlation present in this subsample is still significant (> 99.9%) (see Sect. 3.2). 

Figure 3. The clusters membering the Disk group are plotted in the plane |My | —3.5logRcc 
- logpo • The correlation coefficient is enhanced by ~ 0.2 with respect to the monovariate 
case I My I - logpo ■ As can be seen, the relation is impressively narrowed with respect to 
the monovariate case (Figure 2). This is a further indication that Disk clusters keep memory 
of the "signal" impressed by both \Mv\ and Rqc ■ 

Figure 4. The final cluster distributions in the C -logM plane for the simulations carried 
out in run I (panel a) and run II (panel b). In both cases the imposed correlation {run I) and 
non-correlation {run IT) were preserved till the present epoch (see Sect. 4). The line drawn 
in panel a indicates the relationship between c and logM assumed for the initial conditions. 

Figure 5. The final cluster distributions in the C — log M plane from the run I for inner 
clusters {Rgc < SKpc; panel a) and for outer clusters {Rgc > SKpc; panel b). The line 
drawn in both panels indicates the relationship between C and logM assumed for the ini- 
tial conditions. The correlation appears effectively weakened in the Inner sample and well 
preserved in the Outer one. 
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